The Free Boundary Problem in the 
Optimization of Composite Membranes 

S. Chanillo, D. Grieser, K. Kurata 
February 1, 2008 

Abstract 

In this paper, continuing our earlier article [CGIKO], we study 
qualitative properties of solutions of a certain eigenvalue optimization 
problem. Especially we focus on the study of the free boundary of our 
optimal solutions on general domains. 

1 Introduction and Summary of results 

In this note, we study some qualitative properties of solutions of a certain 
eigenvalue optimization problem. Our note is a continuation and also a 
summarization of the key results of our earlier article [CGIKO]. Our problem 
can be stated in physical terms as : 

Problem(P) Build a body of a prescribed shape out of given materials of 
varying densities, in such a way that the body has a prescribed mass and 
with the property that the fundamental frequency of the resulting membrane 
(with fixed boundary) is lowest possible. 

The physical problem can be re-formulated as a more general mathemat- 
ical problem. More precisely, we are given Q C R n , a bounded domain with 
Lipschitz boundaries and numbers a > 0, A e [0, |0|] (with | • | denoting vol- 
ume). For any measurable set D C Q, let xd be the characteristic function 
and Xa(a, D) the lowest eigenvalue A of the problem, 

— Au + axDU = Xu on Q (1) 
u = on dQ. 
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Define, 

An(a,A)= inf X n (a,D). (2) 

Dctt,\D\=A 

Any minimizer for (§) will be called an optimal configuration for the data 
(ft, a, A). If D is an optimal configuration and u = u a ,D satisfies (0), then 
(u a ,D,D) will be called an optimal pair (or solution). The mathematical 
problem then reads, 

Problem (M) Study existence, uniqueness and qualitative properties of 
optimal pairs. 

We will hereon always work under the nomalization 
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u = 1, u > 0. (3) 

Furthermore, changing D by sets of measure zero does not affect \n(a,D) 
or u, thus sets D that differ by sets of measure zero will be said to be equal. 

A basic tool that we use to analyse our problem is a variational charac- 
terization of eigenvalues, precisely, 

\n{a,D)= mf Rq(u, D), R n (u, D) - 



The minimizer u is well known to exist and is an eigenfunction. Thus, for 
An (a, A) we have 

A n (a,A)= inf Rn{u,D). 

ueH%(n),\D\=A 

The theorem that follows is basic to the questions we hope to treat in this 
paper. The proof of this theorem is to be found in [CGIKO]. To state our 
theorem we will need to introduce some notation. First, we will consistently 
use the notation {u = t} for {x G f2; u(x) = £}, and {u < t} for {x G 
Sl;u(x) < t}. 

Theorem 1 ([CGIKO]) For any a > and A G [0, there exists an 
optimal pair. Moreover any optimal pair (u a> D, D) has the following proper- 
ties: 

(a) Ua,D G C 1 ' 5 ^) n W 2 ' 2 {Vt) n C 7 (H) for every 5 < 1 and some 7 > 0. 

(b) D is a sub-level set ofu a ,D, i-e. there exists t > such that 

D = {x G Q; u a) n(x) < t}. 
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(c) The value of a for which An (a, A) = a, is unique. 

(d) Every level set {u a ^o = s},s ^ t has measure zero. If in addtion 
An(a, A) ^ a, the free boundary, T = {w a ,D — t} has measure zero. 

We use the notation Uq(A) for the unique value of a in part (c) above, 
that is, 

a n (A) = A n (an (4) 

We also see right away that our problem to determine an optimal pair 
(u a ,D, D) which seemed linear is now a non-linear problem, 

— A.u a ^D + Qzx{u atD <t} u a,D = A n (a,A)u a , D on Q (5) 

u a ,D = on dfl. 

Another important remark is that because any optimal configuration D is a 
sub-level set and u a ^ = on dQ, the set D will always contain a tubular 
neighborhood of dQ, i.e. D always contains a boundary layer. 

For notational convenience, from now on we will drop the subscript Q 
and write A(a,A) for An(a, A), and use the simpler notation u for -u Q) d, in 
situations where no confusion arises. || • ||oo will denote the supremum norm 
in L°°(Vt) and || • H2 the norm in L 2 (Q). 

The main focus in this paper will be on the free boundary {u = t} on a 
general domain Q. 

Before we state the theorems that we prove in this paper, we continue 
summarizing some of the salient results of [CGIKO]. It is proved there that 
problem (M) generalizes problem (P). In fact for a < a^(A), the solutions 
of problem (M) and (P) are in one to one correspondence. Another natural 
questions that arises is if D inherits natural symmetries that fl possesses, and 
if given Q, does there exist a unique optimal configuration D. The answer is 
negative on general domains unless Q has very strong topological restrictions. 
Symmetrization and rearrangement invariant integral methods allow one to 
prove the next theorem. 

Theorem 2 ([CGIKO]) Assume Q is symmetric and convex with respect 
to the hyperplane {x x = 0}. That is for each fixed x' = (x 2 , ■ ■ ■ , x n ) the set 

{xi : (x 1 ,x') E f2} 

is either empty or an interval of the form (— c, c). Then any optimal solu- 
tion (u, D) is symmetric with respect to the hyperplane {x\ = 0} and u is 
decreasing in x\ for xi > 0. 
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Theorem |2] implies the next corollary, which is the only uniqueness result 
proved in [CGIKO]. 

Corollary 1 Let Q = {\x\ < 1} be the ball. Then the optimal configuration 
is unique for any a, A and furthermore D is an annular region, 

D = {x; r(A) < \x\ < 1}. 

On other domains we encounter the phenomenon of symmetry breaking. 
Specifically in [CGIKO] we show symmetry breaking phenomena on annular 
domains in R 2 and on dumbbell shaped domains. We have 

Theorem 3 ([CGIKO]) Fix any a > and 5 G (0,1). Let Vl a = {x e 
R 2 ; a < \x\ < a + 1}. Then there exists ao = a (a,5), such that whenever 
a > ao and D is an optimal configuration for Q a with parameters a and 
A = 5\Q a \, then D is not rotationally symmetric. 

Because Q a is rotationally invariant, Theorem ||] implies that there are in- 
finitely many choices for the optimal configuration D on annuli. On dumb- 
bell shaped domains we have, in addition to symmetry breaking, some extra 
information on D. We define the dumbbell shaped domain J\ by 

n h = £?i((-2, 0)) U B 1 {{2, 0)) U ((-2, 2) x (-h, h)), (6) 

where B r (p) = {x E R 2 ; \x — p\ < r}. We call the disks B r (p), the lobes of 
the dumbbell and the strip (—2, 2) x (—h, h) the handle. We have 

Theorem 4 ([CGIKO]) For any given a > and A e (0, 2ir), there exists 
h > such that for domains £\ of (^) with h < h , 

(a) Any optimal pair (u,D) is not symmetric with respect to the x^ — axis. 

(b) If A > n, then for any optimal pair (u,D), D c is totally contained in one 
of the lobes ^((±2,0)). 

We end our summary of results from [CGIKO] with a theorem on convex 
domains. 

Theorem 5 ([CGIKO]) Suppose Q is convex and has a smooth boundary. 
Then there exists ao(Q,A) > 0, such that for any a < oq and any optimal 
configuration D, one has 

(a) dD n f2 is real- analytic. 

(b) D c is convex. 
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Theorem [| should be compared with the basic theorem of Brascamp-Lieb 
[BL] that establishes the convexity of level lines of the first eigenfunction on 
convex domains, when a = 0. Theorem |5| extends the result of [BL] to some 
values of a > 0, but it is completely open if Theorem || extends to the case 
of all a > 0. 

We now turn to the results that we prove in this paper. Our focus will 
primarily be on general domains and in particular on the free boundary for 
the optimal pair (u,D). In a general domain Q the first eigenfunction if) 
(with standard L 2 normalization) that is, 

— Aip = fixip on Q (7) 

ip = o on on, I if? = 1, 

is real-analytic in the interior and this places very strong restrictions on the 
exceptional sets, i.e. places on the level sets of if), where = 0. If we were 
in R 2 , the exceptional set would consist of points. In this analysis unique 
continuation plays a role, since we easily see that if w = ip Xi then from (0), 
—Aw = fi\W, and thus unique continuation yields some information on the 
zero set of w. Thus in our problem it is clear the free boundary {u = t} 
on a general domain will possess an exceptional set, since if) in general has 
one, but any attempt to understand the fine structure of the exceptional 
set, Hausdorff measure, rectifiability etc., through a unique continuation ap- 
proach is difficult. The reason being, unique continuation will not apply, 
since u is only weakly regular. This prevents us from obtaining an equation 
for W\ = u Xi . In addition a further difficulty is that (|5|) is an equation of 
the type — Au + V(x)u = 0, and this prevents us from obtaining a homoge- 
neous equation that is satisfied by w 2 = u — t to which we may apply unique 
continuation to study the level surface {u — t},t^ 0. 

Another approach is to view our problem (|5|) as a perturbation in a from 
the problem (|7|). This approach suffers from the fact that we do not get 
additional information for large a. One may view this again as a difficulty 
arising from lack of continuation properties in a. Two results in this direction 
are proved in [CGIKO]. We reproduce the statements and proof here. 

Theorem 6 For s > 0, let [Q] s = {ip < s}, where ip is the normalized 
first eigenfunction of problem (Q). Fix A G [0, and choose tn such that 
|[f2]' n | = A. Then for any 5 > 0, there is ao = ao(8, fl) such that if a < cto 
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and D is an optimal configuration for (a, A), then \t — tn\ < S and 

pyn-s cDc [Qf»+ S . 

The basic lemma that is used to prove Theorem |6| can be used to analyse the 
limiting behavior as A — > |0|. We have 

Theorem 7 Let Q be a smooth bounded domain. Let a > be fixed. Let ip 
be the function of (Qj and let 

M = m&xib. 
n 

Then for any 5 > 0, there is A Q = A Q (5, a,Q) < \Q\, such that whenever 
A > A and D is an optimal configuration for (a, A), then 

D c c {ip > M — 5}. 

The meaning of Theorem |6| is that the free boundary for our optimization 
problem, that is the set {u = t}, is "trapped" between the levels t^ — 5 and 
tn+5 of the first eigenfunction ip for the domain Q. However, this information 
is too weak to conclude anything fine about the free boundary even for small 
a > 0. Theorem [F]on the other hand indicates that as A — > |0|, D c coalesces 
onto the set where ip achieves its maximum, ip being the first eigenfunction 
on Q, see (^). Now, keeping in mind part (b) of Theorem [|, it is likely that D c 
may coalesce onto a strict subset of the set where tp achieves its maximum. 

Even though we have been unable to apply unique continuation to study 
the free boundary for large a, it is still possible to apply the Hopf lemma 
[GT, Lemma 3.4] and get some information on the free boundary. 

A typical result we prove is: 

Theorem 8 Let a > A(a,A). Let T = {u = t} denote the free boundary 
set. Then there is a subset £ of T such that 

(a) £ is a G$ set. 

(b) T\£ is a real- analytic, n — 1 dimensional sub-manifold ofH n . 

(c) If moreover a > A(a,A), then for every xq G T and every e > 0, the 
ball B t (x ) contains points of both {u > t} and {u < t}. 

We refer to the set £ as the exceptional set. From the construction of the 
exceptional set, we will deduce further geometric information regarding the 
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free boundary. This is the content of Proposition 0, which we do not state 
here (see section 2). 

We now give a sufficient condition that ensures that the hypothesis of 
Theorem |a> A (a, A) is fulfilled. 

Proposition 1 Let a > /ii(fi), where is the first Dirichlet eigenvalue 

for —A on Q. Then there exists Aq = Ao(a), such that a > A(a,A) for all 
A < Aq. Furthermore, for Ci = \\ip\\^ and for fixed A e (0, C\), there exists 
ao such that a > A(a, A) for all a > Oq. 

We lastly investigate the effect of the curvature of dQ on the free boundary 
and the "thickness" of the optimal configuration. As observed earlier in the 
remarks after Theorem [TJ, D always contains a tubular neighborhood of the 
boundary. The theorem that follows demonstrates in the model case of an 
annulus in R 2 , that at places where dfl has large "negative" curvature, one 
finds that D is "thin" . To state our result we will set up some notation. Let 

Q t = {x e R 2 ; e < \x\ < f }, B = {x G R 2 ; \x\ < 1}. 

For fixed e , let A < 7r(l — e^). For e < e and any fixed a > 0, let (u e , D e ) 
denote the optimal pair for Q e , with lowest eigenvalue An e (a,A), with con- 
straint \D e \ = A. We now claim that for ni(B), the first Dirichlet eigenvalue 
of the unit disk, we have 

An e (a,A)>//i(S). (8) 

Since Q e C B, using (u e , D e ) as a trial pair in the variational characterization 
for Aci e (a, A), we have 

An e {a,A)= [ \Vu e \ 2 + a [ \bA> [ |V^| 2 > fMi(B). 
Jn E Jn e jo £ 

This establishes, (||). As a consequence of @, imposing the hypothesis a < 
fJ>i{B), ensures that for every e > 0, a < an e (A) (see the definition (||) ). 
Next, if a < aa € (A), and if D e is radially distributed, Theorem 2 of [CGIKO], 
yields that D t has the form, 

D e = {x E R 2 ; e < \x\ < r € or R € < \x\ < 1} (9) 

for some r e , R e , e < r e < R e < 1. Thus if a < fi\(B) we may assume that if 
D € is radially distributed, then the set D e has the form described by @ for 
every e > 0. We have 
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Theorem 9 Assume a < fi\(B), and A > is prescribed. Given this choice 
of a and A, let (u e ,D e ) be an optimal pair with \D e \ = A. Assume D t is 
radially distributed and hence of the form (j^j. Then, 

limsupr,, = 0. 

e->0 

Thus the implication is that D e thins out on the boundary layer in contact 
with {x; \x\ = e}, the inner boundary of <9f2 e . As e — > 0, the curvature of 
{x; \x\ = e} is increasing and negative as seen from Q e . Thus in this model 
case one may conclude that diam (D) is small on parts of D which are in 
contact with pieces of dfl, where the curvature of dQ is large and where dfl 
is concave as seen from Q. 

The paper [CGIKO] discusses the historical antecedents of this problem 
and the interested reader is referred to this paper for a discussion. Further- 
more the optimization problem discussed here, is amenable to being modelled 
on a computer. Details of the numerical simulation are available in [CGIKO] 
and the interested reader may find the source of the algorithms used and the 
shape of the optimal configuration in many types of domains obtained by 
these numerical studies. 

2 Proofs of the Theorems 

In this section we prove Theorems |6|-|9] and Proposition |l| We begin with 
the proof of Theorem ^|. We need a preparatory Lemma, that is well-known 
in perturbation theory and in the Physics literature [B, Appendix 39, p. 
469] . We want a slightly more precise form, though the technique of proof is 
standard and the basic idea follows from [B]. 

Lemma 1 Fix D C Q. Let u a ^ be the ( positive, L? normalized ) first 
eigenfunction of —A + axD with eigenvalue A (a, D). Then there is a constant 
C = Cq such that for < a < 1 (ip, fa refers to ^)), 

(a) < A(a, D) — fa < a, 

(b) \\u a>D - 4>\\Hi(n) < Ca, 

(c) \\Ua,D ~ lf>\\L°°(Cl) < Ca. 

Proof of Lemma 1: Recall we have set u a ^ = u. Note 

X(a,D) < /(|V^| 2 + axD^ 2 ) < fa + a(^ 2 
Jn Jn 



8 



< fi\ + a. 

Thus, A(a, D) — fix < a. Next, 

\(a, D) = / (| Vw| 2 + axDU 2 ) > \i\ + a / Xdu 2 . 
Jn Jn 

Thus, X(a, D) — \i\ > 0, and we have (a). To prove (b), let {ip^kLi be an 
orthogonal basis of eigenfunction of —A with Dirichlet boundary conditions 
(Note ipi = ip of problem ([])). The corresponding eigenvalues will be denoted 
by {iik)T=ii where it is well-known that \i\ is simple. Expanding u, we have 
u — Y^jLi Pji/ij , and thus (—A — jii)u = Z)J1 2 /%(/•*:? ~ ^1)^3 an< l 

IKA + ^II^E^^-^) 2 - (10) 

where || • H2 denotes the L 2 (f2) norm. From — Au + cuxdu = A(a, D)u, we get 

(—A — ni)u = (A(a, D) — ni)u — oixdu. 

Therefore applying (a), ||(A + /ii)u||2 < Ca. Since \i\ is simple, there exists 
5 > 0, 5 = 5(Q), such that fij — fii > 5 > for j > 2. Then from fllTf ) we get 

Ca 2 

I 2 



E^<^- (id 



i=2 



We re-write u as u = Piipi + and (|TT]) gives ||^||2 < Cc^ X - Now 1 
Hh = Pi + thus 



l + /3i — 5 2 

Here we used the fact that /3i = («, ^l) > 0, because both u and t/>i are 
positive. Therefore, 

Co? 

\\u ~ -0x111 = (A - I) 2 + < ^- < Ca\ (12) 
All the remaining consequences follow from flT2"|). From (|5[), 

- A(u - ^1) = (A(a, .D) - yUi)w + yUi(w - - oxdM- (13) 
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We re-write (|T3|) as 

— A(u — ipi) — (jli(u— ipi) = (A(ot, D) — fMxju — oxdu = g. 
Now from [GT, Theorem 8.15] again, it follows that 

\\u - ^llloo < C\\u - Vi II 2 + Ca- 
using ([12]) on the right side, 

\\u - ^i||oo < Ca, 

which is part (c). Using ||g||oo < Ca and part (c), we conclude ||A(u — 
^i)||oo < Ca, which is (b). □ 

Theorem |6| and [7| are now consequences of Lemma |]. 
Proof of Theorem |6|: Apply Lemma [I] (c) to the optimal pair (u,D). 
Choose a = 5/{2C), so that \\u — ipi\\oo < 8/2 f° r a < «o- From Theorem [1], 
if x G D,u(x) < t, and so ip(x) < t + 5/2 and hence D C [f2]* +<5/2 . In a 
similar way we establish [fi]' -<5 / 2 C D. Thus we have 

pt-«/2 c £, c [O]^ 5 / 2 . 

From the statement above, we get |[fi]' _<5 / 2 | < A < \[Q] t+s / 2 \, and thus by 
continuity, there exists ta such that A = \[Q] tn \, and \tn — t\ < 5/2. From 
this assertion the assertions of Theorem ^ follow. □ 

Proof of Theorem [7|: We begin by showing that a slight modification of 
the proof of Lemma [I] yields, 

||tt-Vi||oo<C r a , n (|fi|-A). (14) 

We show first, 

- (A(a, A) - a)| < Ca,n(|n| - A). (15) 
We re-write our equation for u as 

— Au — axD^u = (A — a)u, A = A(a,A). (16) 

From (|l^) we have 

Hi — a Xd^u 2 < / | Vw| 2 — a / Xd^u 2 = A — a. 
Jn Jn Jn 
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Thus, 

//i - (A - a) < a [ xdcu 2 < C\D C \ = C(\Q\ - A). 
Jn 

Next, we have 

A - a < (|V^| 2 - axDcip 2 ) < fix -at / Xd^ 2 
Jn Jn 

which yields < \L\ — (A — a). The assertion (|15|) follows. Using (|15|) 
and the equation (|i~6l) we can proceed as in Lemma [l| to obtain (|14|) . If 
A< 6/(2C a>n ), from fllj) we see [ft]*- 5 / 2 C D = {u < t} C [ft]'+ 5 / 2 . 
Since \D\ = A, we have A < |[ft]' +5 / 2 | = \{ij> <t + 5/2}\. Thus if in addition 
A > Aq, we can arrange the situation so as to have M — t < 8/2. So 
{tt] M - s C [ny- 5 / 2 C .D. The conclusion [n} M ~ 5 C D is readily seen to be 
equivalent to the assertion made in Theorem [7| □ 

We need some preparatory lemmas before we prove Theorem |8|. As usual 
u = u a,D w iU denote the solution to our optimization problem and conse- 
quently u will satisfy (H). 

Lemma 2 (a) Fix any a > 0. Let the free boundary set be T , T — {u — 
t}. We let D + = {x;u(x) > t}. Assume D + satisfies an interior sphere 
condition with respect to x G T , that is there exists a ball B, B C D + and 
0BnJ r ={x }. Then\Vu(x )\ ^0. 

(b) Let D~ = {x; u(x) <t} (in the situation a ^ a^{A), D~ = D). Assume 
that a > Aq(oj, A). Let D~ satisfy an interior sphere condition with respect 
to xq G T , that is there exists a ball B, B C D~ and dB n T = {xq}. Then 
\Vu(x )\^0. 

Proof of Lemma The proof of both parts of our lemma rely on Hopf 's 
lemma [GT, Lemma 3.4]. We prove (a). Set <ft = t — u. We observe that in 
the ball B C D + , u satisfies from (^) 

-Am = A(a, A)u. 

Thus A0 = A(a, A)u > in B, and < on B with 0(x o ) = 0. Hopf 's 
lemma then yields |V</>(xo)| = |Vm(xo)| ^ 0. 

(b). The proof of this part is similar to part (a). Since B C D~, from 
(H) we see on B we have 

— Au + au = Au. 
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Since a > A, we easily see Au > on B. Thus on B G D we have 
Acf> < 0, > on B and 4>{xq) = 0. Thus Hopf's lemma again yields 
|V0(x o )| = \Vu(x )\^0. □ 

Lemma 3 Let h(r],p),r] G R, p = (pi, • • • ,p n ) £ R™ ^ e ° locally bounded 
function. Let w G C 1 ^) satisfy 

Aw = h(w,Vw). (17) 

Assume furthermore h is smooth in the variable p. Assume at the point 
x G Q, Vro(xo) 7^ 0. T/ien there exists a ball B, x G B, such that the set, 

{x G B; w(x) = w(x )} = S 

is a smooth hypersurface ofH n . If in addition h is real-analytic in the variable 
p, the set S is also real- analytic. 

Proof of Lemma 3: Since h{w, Vu>) is locally bounded, it follows by elliptic 
estimates that w G C 1 ' 7 fl W 2 ' s , s < oo. Thus by the implicit function 
theorem, since Vw(io) ^ 0, we conclude that S is a C 1 ' 7 hypersurface for 
all 7 < 1. Now we shall improve the regularity of the hypersurface S. By 
a rotation of coordinates we may assume w Xi (xq) = 0,i = 1, • • • ,n — 1 and 
u! Xn (xo) 7^ 0. Let x' — (x±, • • • , a? n -i) and consider the map, 

^ : 5 -> R", V(x', X n ) = (X', W(X', X n )). 

We denote points in the image of by y = (y', y n ) where y' = (yi, • • • , y n ~\) 
and y n = w(x', x n ). Let ^Z^ 1 denote the inverse map to which will exist if 
B is picked to be small. We have, 

*-\y',y n ) = {y',F{y',y n )). 

Now, 

F(x',w(x',x n )) = x n . 
Differentiating the equation above we get the equations, 

Fy.+w Xi F Vn = 0,i = l,---,n-l, and F Vn w Xn = 1. (18) 
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By the chain rule, 



d d d d d 

+ w Xi - — , - — = w x 



dxi dyi l dy n dx n "dy n 
From (0), w Xi = -F y jF yn , i = 1, - • • , n - 1 and w Xn = l/F Vn . Thus, 

d d \ f-F v .\ d ( 1 

£i v% dy n J \ F yn J dy n \F Vn 

n ^(± _ ZkJM ( z*k \ + J_ JL( J_ 
^ V % F,„ c>y n )\F y J F yn dy n \ F yn 

= LF. 

Next we freeze the coefficients of L at xq. Since w Xi (xo) — 0, i — 1, ■ ■ ■ , n — 1 
and w Xn (^o) = a 7^ 0, we see from (|18|), F y4 = and i 7 ^ = a -1 at \I/(a:o)- At 
yo = \I/(xo), we have 

1 ri ~ 1 1 
rp = _ Vp —F 

±J± t? 2/»2/i p3 



F, 



n— 1 j 



UnUn 



Vn L i = l - 1 y 



Thus if S is picked suitably small, LF is an elliptic, quasi-linear operator. 
Since w satisfies (|17D, F satisfies 

LF = h(y n ,A{=j^,^-)) (19) 

Vn Vn 

where A is a fixed matrix in 0(n), the rotation group, associated with our ro- 
tation of coordinates, and A(zi, z n ), denotes the product of A and the vector 
(zi, ■ • • , z n ), (zi) = (zi, • • • , z n -i). Since w G C 1 ' 7 fl W 2 ' s for all 7 < 1, s < 00, 
we see the coefficients of L are in C ' 7 fl W 1,s , s < 00. Differentiating 
equation ( |T9"D in any of the variables y iy i 7^ n, we may apply a standard 
bootstrap argument using well-known elliptic estimates for example [GT, 
Theorem 9.11] to conclude from the fact that h is smooth in the variable p, 
that F is smooth in any of the variables yi, i — 1, • • • , n — 1. If one has in 
addition that h is real-analytic in the variables p, then applying the results 
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of Morrey [M, Theorem C] or Friedman [F, Theorems 1, 4], we can conclude 
that F is real-analytic in the variables yi,i — 1, ■ ■ ■ ,n — 1. Now the defining 
equation for S is given by x n = F(x\, ■ ■ ■ , x n _i, w(x )). Since F is smooth 
(real- analytic) depending on the regularity of h in the p variables, it follows 
that S is a smooth (real-analytic) manifold, depending on the fact that h is 
smooth (real-analytic) in the p variables. □ 
We are now ready to prove Theorem ||. 
Proof of Theorem ^ We construct the exceptional set £. Let 

K n = {i G fi; distance (x,J-) = — }. 

n 

Now define ^ 

F n = {x G T\ distance (x,K n ) = —}. 

n 

The sets K n and F n are closed sets for all n G N. We define the exceptional 
set by 

e = r\{\j^ x F n ). 

Since each set T fl F£ is open in J 7 , £ is a Gs set. This proves Theorem || 
(a). 

We shall now prove that for each point Xq G J-\£, we can construct 
either an interior ball B C Z^ + , such that dB fl JF = {xo} or an interior ball 
B C -D - , such that dB fl JF = {x }. Since we are assuming a > A(a,A), 
Lemma ensures that for each x G T \ £, Vu(xq) ^ 0. Now the PDE 
satisfied by u, that is (|5]), can be written as 

Au = h{u), 

where h{rj) = —A(a,A)r] + axG( 7 l) r l^ = {ViV < t}- Thus h G L^ C (R) 
and the hypotheses of Lemma [3] apply to u. Thus there exists a ball B , 
centered at Xq, such that, T fl So is a real analytic manifold. So we now 
verify our claims regarding the interior spheres. Fix a point xq G T \ £. 
Then x £ F n for some n. Let G -ft'n, such that \xq — z \ = 1/n. We claim 
the ball Bi/ n (z ) is totally contained in D + or D~ . Suppose there are points 
Z\ G D + and z 2 G -D~ in B\/ n (zo). Then by the continuity of u, the line 
segment joining z\ to ^2 which also lies inside B\/ n (zo) will contain a point 
of J 7 . Thus distance (zq,^) < 1/n and hence zq G" -R'n- Now pick a ball 
5', 5' C By n (zo), and centered along the radius joining xq to zq and with 
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dB' D dBy n (zo) = {xo}- Clearly B' \ {x } is contained either in D + or D~ 
and dB' n T = {xq}. The hypotheses of Lemma ^| are fulfilled. Theorem ||] 
(b) now follows. 

To prove Theorem |^ (c), we can argue via the strong maximum principle. 
Since a > A(a,A), applying Theorem [I] (d), we see that D = {u < t}. If 
for some e > 0, the ball B e (xo), Xq G T contains no points of D~ and only 
points {u > t}, then on B e (x ), u satisfies 

-Am = A(a, A)u > 0. 

So Au < on B e (xo) and u > t on B e (xo) with u(xq) = t which is a 
contradiction. We may argue as in Lemma |2], part (b) and show on B e (xo) 
there are also points of D + for every e > 0. □ 

To discuss further geometric properties of S, we introduce, 

K e = {x G f2; distance (x, T) = e}, 

F e = {x G T\ distance (x, K t ) = e}. 

We have 

Proposition 2 For every a > 0, 

(a) F ei C F e2 fort! > e 2 . 

As a consequence of (a), we have 

(b) U n >iF n = U €>0 F e , F n C F m for m> n, andE = T\ U e>0 F e . 

(c) If a > A(a,A), then for every point x G £, any ball B e (x ) contains 
points y + G D + and y_ G D~ such that distance (y±,J r ) < \y± — x \. 

(d) Furthermore, if z G T , such that \y + — z \ = distance (y + ,J r ), then D + 
satisfies an interior sphere condition in the sense of Lemma |^ with respect 
to Zq. Thus by the proof of Theorem there is a neighborhood B of zo, 
such that {x G B : u(x) = u(z )} is a real-analytic hypersurface. A similar 
statement holds for y_. 

The meaning of (b) is that since £ = T fl (fl^L 1 F^) and F£ D F^ for m > n, 
the exceptional set is really a G$ set formed by the intersection of the nested 
open sets F£. The meaning of (c), (d) is that the behaviour of the free 
boundary in the neighborhood of the exceptional set at least in R 2 is that of 
isolated singularities with a conical structure, with the cone having its vertex 
at xo G £ . The cone locally divides R 2 into at least two components, one 
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component is contained in D + and the other is contained in D~ . This geo- 
metric picture is only heuristic since it still needs to be rigorously established 
that the component of T that contains Zq also contains Xo G £. Only then 
can we conclude there is a true conical singularity at Xq. 

Proof of Proposition |2|: Fix x G F ei . Then by definition, one can find 
z G K ei , such that \z — x\ = e\. Since z G K €1 , the ball B €1 (z) will contain 
no points of T in the interior. Now consider the radial line connecting z 
and x, and locate on this line a point y such that \y — x\ = t^- The ball 
B €2 (y) C B ei (z), and x is the sole point in T on dB e2 (y). Now by definition 
y G K t2 , and since x G T and \x — y\ = e 2 , x G F e2 . This proves (a), and (b) 
is then an elementary consequence of set theory. 

Next, by Theorem |3| (c) we know that the ball B € (xq) with xq G S ( in 
general for any point in T actually), contains points y± in D ± respectively. 
Let \y + — Xq\ = 5. We claim distance (y + ,J r ) < 5. If distance (y + ,J r ) > 5, 
then y + G K T , r = distance (y + , J? 7 ), r > 5. Since \y + — x \ = 5 < r, it follows 
that x G F T . Thus Xq <^L S. This proves part (c), since an argument similar 
to the one above takes care of the point ?/_ G D~ . 

Lastly we prove part (d). Now, \y + — zq\ = r, zq G T . Thus by the argu- 
ment employed in Theorem |8|(b) it is easily seen, that the open ball B T (y + ) 
contains only points of D + . Again employing the argument of Theorem |8] 
(b), we can find a ball B' C B T (y + ) such that dB' n T = {zq} and hence B' 
is the desired interior ball. □ 

We now prove Proposition [[]. 
Proof of Proposition [l]: Arguing as in Lemma 0(a), but now using the 
fact that Halloo < C, [GT, Theorem 8.15], we have 

A(a,A)< f \V^\ 2 + af X d^ 2 < fii + aW\lA. 

Thus, 

A(a, A) < + aC A, C„ = |hC (20) 

Since a > we can find A > such that aCoA Q < a — It follows from 
( |20D that for A < A Q , a > A (a, A). The second part of Proposition [I] also 
follows from ©. Select C x = C _1 . If A < C h C A = 1 - e, e > 0. Thus for 
a > y^i < and hence by (|20|), A(a, A) < a. □ 

Proof of Theorem ||: We set A e = A Uc (a, A), A = A B (a, A). D will denote 
the optimal configuration for B, that is D = {x; r < \x\ < 1}, 7r(l — Tq) = A. 
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We also need to consider the first Dirichlet eigenvalue X e , of the problem, 

-Ag + axDg = Kg in Q e 
g\an e = 0. 

We claim 

0< A e -A < Clloger 1 . (21) 

Extending u € to {|x| < e} by setting u t = on {\x\ < e} and using this 
extended function trial function with D f cLS cL trial configuration on 

B = {\x\ < 1}, we see A < A e and so we have the left side in (21). Next by 
Theorem 2 in Swanson [S], 

< A e - A < c| logep 1 . (22) 

In fact < A e — A is simply a consequence of domain monotonicity. By the 
variational characterization A e < A e , and thus from ( 22f) we easily have (|21|). 

We now establish limsup e ^ r e = 0, by contradiction. Assume there is a 
sequence €j \ 0, lim^oo r e . = 5 > 0. Then R e — > b$ and since \D e .\ = A, 
the limit set T>s is: 

Us = {x; \x\ < 5 or b$ < \x\ < 1}, \U$\ = A. 

We set, 

D$,e = {x; e < \x\ < 5 or b$ < \x\ < 1}. 

We use the notation X e (Ds >e ) for the first Dirichlet eigenvalue on Q e for the 
problem, 

-Aw + axD^w = X e (D Sj e)w on Q e (23) 
w = on dQ e , / w 2 = 1. 

X(Us) will denote the first Dirichlet eigenvalue for —A + a\v s on B. We 
claim, 

|A e . -A(P*)|->0 as e.-^O. (24) 

We have 

\K e .-X{U s )\ < |A ei -A ei (^)| + |A e3 .(^, ei )-A(P 5 )| 
= Ji + J 2 - 
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By using Theorem 2 in [S], one can conclude J2 < C\ logej| _1 — > as — >■ 0. 
We now show 3\ — > 0. Let (w e , D e ) be the optimal pair for fl e . Then, using 
the eigenfunction w from (pf), and the uniform bounds ||w||.L°°(n 6 ) < C, C 
independent of e > 0, which follows from [GT, Theorem 8.15], we have 



A, < 



/ {\Vw\ 2 + a XDt w 2 ) 

= / (\Vw\ 2 + axD 5 ,w 2 ) + a (XDc-XD s ,Jw 2 

< X e {D s>e )+Ca\D e AD S!e \, 

where D e A D$ e is the symmetric difference of the sets D e , D§ € . Thus for the 
sequence 6j, we easily have \D e A Ds >e \ — > as j — > 00. We conclude 

A ej <A e . (Di, ej )+o(l), j — > 00. 

Likewise using the function u e in the argument above, we have 

K s (D s>ej ) < A ej +Ca\D €j AD s>€j \ 
< A e . +o(l), j — > 00. 

Thus, 

\X ej (D 6>e .) - A £j \ ^ as j -»• 00. 

Thus Ji — > as €j — > 0. This establishes fl24l) . We infer from ([H]) and (0) 
that as j — > 00, 

I A - A(Z>,)| < I A - A e .| + |A e . - \{V S )\ - 0. 

Hence A = \(T>s). However, this contradicts Corollary [l] of our introduction, 
a proof of which is supplied in [CGIKO] . □ 



3 Open Problems and Conjectures 

A number of open problems and conjectures can be stated based on the 
numerical data in [CGIKO] and the rigorous results there and also on results 
proved here. We will outline some. 

Problem 1 : (Uniqueness of the optimal configuration) The only domain 
for which we have established the uniqueness of the optimal configuration is 
the ball, see Corollary Is D unique if Q is convex? 
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Problem 2 : (Continuation Problem) Theorem^ states that on a convex do- 
main Q, D c is convex for small a > 0. Is it possible to continue along the 
values a > 0, to obtain convexity of D c for all a > 0? 

Problem 3 :(The free boundary on general domains) The free boundary T 
on general domains will contain an exceptional set £ as constructed in the 
proof of Theorem 8. What is the Hausdorff dimension of S? Is at least T 
a rectifiable set? One suspects S consists of points, where real-analytic arcs 
intersect if Q C R 2 . 

Problem 4 : (Monotonicity of D) Suppose A < A', then does this imply 
D a A C D a A> ? If symmetry breaking occurs this statement needs to be mod- 
ified. Nevertheless on domains where the optimal configuration is unique, 
does the above monotonicity hold? 

Problem 5 : (Symmetry breaking on annuli) When Q is an annulus, what 
is the shape of D precisely? The proof of Theorem^ in our introduction as 
presented in [CGIKO] and the numerical computations presented in [CGIKO] 
suggest that D c lies between two rays 6 = and 9 = (3. In fact the results 
in [CGIKO] suggest that, (3 = n/N,N = N(a,\D\/\Q\), and N -»• oo as 

\d\ -> \n\. 

Problem 6 : (Influence of the boundary curvature) In Theorem we saw 
in a model case that the diameter of D is affected by the curvature of dVt. 
Investigate this phenomena on general domains fl. 

We refer the interested reader to [CGIKO] for further problems and conjec- 
tures. 
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